Pore Stability and Dynamics in Polymer Membranes 
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Vesicles self-assembled from amphiphilic diblock copolymers exhibit a wide diversity of behavior 
upon electroporation, due to competitions between edge, surface and bending energies that drive the 
system, while different viscous dissipation mechanisms determine the timescales. These effects are 
manifestations of the varying membrane thickness d from what are essentially chemically identical 
molecules. For smaller d, we find large unstable pores and the resulting membrane fragments 
reassemble into vesicles within minutes. Vesicles with larger d form long-lived pores of nanometer 
dimension that could be potentially useful for sieving applications and drug delivery. 

PACS numbers: 82.70.Uv, 68.65.-k, 87.68. +z 



Membranes are ubiquitous in biology and play a central 
role in the distribution of macromolecules both inside and 
outside cells. Additionally, specialized transmembrane 
machinery serves to exquisitely regulate intracellular con- 
centrations at their optimal, non-equilibrium conditions 
jjj. These transmembrane gradients are known to be in- 
volved in signaling and other regulatory processes. Thus 
the control and exploitation of biomembrane selectivity 
has been and continues to be the focus of intense research 
Q . One approach to manipulate biomembranes is known 
as electroporation, where an electric field is used to tran- 
siently disrupt the membrane ||. In the brief time the 
membrane is open, solutes or macromolecules can diffuse 
into or out of the cell. Technological applications for 
DNA delivery and creation of hybrid cells have employed 
electroporation, although cell viability is not always pre- 
served ||. 

Necessarily generated during biological processes are 
transmembrane potentials V m and membrane tensions r, 
as a result of electrostatic or mechanical forces. The con- 
tinuing obstacle to approaching biological systems is that 
their complexity can obscure underlying physical mech- 
anisms. Taking a reductionist approach, we have cho- 
sen to use synthetic vesicles formed from diblock copoly- 
mers, called polymersomes. These self-assembling di- 
blocks are amphiphilic, with the hydrophilic portion be- 
ing poly(ethylene oxide) (PEO) and the hydrophobic part 
composed of either polybutadiene (PBD) or its saturated 
form, poly(ethylethylene) (PEE) {§. The synthetic ba- 
sis allows for creation of a family of varying molecular 
weight M„ from w 4 — 20 kg/mol [B. The resulting hy- 
drophobic thickness d ranges from 8 — 21 nm, compared 
to biomembranes of 3 — 5 nm. Systematic investigations 
of electroporation are therefore motivated from both a 
biological and a physical perspective. In this Letter, we 
report our observations of a wide diversity of pore behav- 
ior, owing to the expanded range of d that is available 
with polymeric systems. 

The experimental procedure closely follows Aranda- 
Espinoza, et al. In all cases, we used the film re- 



hydration method to prepare vesicles, as is common for 
liposomes. A combination of phase contrast and fluo- 
rescence microscopy were used to monitor any changes 
in vesicle integrity that occured following pore formation 
§• 

The energetics of a circular hole of size r in a flat, in- 
finitely thin, membrane can be most simply described by 
E = 2rrrT — 7rr 2 E, where T and £ are the line and sur- 
face tensions of the membrane, respectively ||. The line 
tension accounts for the energetic penalty involved with 
having a hole and the surface tension reflects the energy 
associated with a loss of membrane area. This relation 
predicts a metastable pore size r* = T/T, obtained from 
dE/dr = 0. Pores of size r < r* will reseal and those 
with r > r* will grow without bound. 

More complex models of pore growth include a dy- 
namic surface tension, predicting a stable pore size 
, , but our observations indicate unstable pores in al- 
most every case for OE7 (d = 8 nm) and OB16 (d = 11 
nm) , indicating a small value of V relative to the surface 
tension S. However, we need to consider the influence of 
membrane thickness on both T and E. At the moment of 
poration, £ is the lysis tension £ c , and previous work has 
shown that E c ~ d |(|. The dependence of T on d, con- 
sidering both hydrophobic and hydrophilic pores, leads 
to the same scaling of T ~ d. For a hydrophobic pore, 
the line tension is most simply the product of the exposed 
length d and the interfacial energy density 7: T ~ jd. For 
a hydrophilic pore, T should scale as the bending modu- 
lus k c multiplied by the curavture ~ 2/d. Simple elastic 
models predict k c ~ d 2 [jy] and hence r ~ d again. 

We therefore have the counterintuitive result that the 
metastable pore size does not depend on membrane thick- 
ness. We do indeed see large pores from OE7 to OB16 
but OB18 (d = 15 nm) and OB19 (d = 21 nm) exhibit 
significantly smaller, and long-lived, pores. This abrupt 
change is presumed to be a result of increasing viscosity 
or chain entanglements within the hydrophobic region of 
the membrane j|, |l2). The lifetime of pores is largely dic- 
tated by the dynamics of the membrane. Liposomes are 
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known to reseal quickly, with a few notable exceptions 
|H . Under no external stresses, viscous dissipation at 
the interface and within the membrane are the obstacles 
to in-plane material rearrangement. 




Figure 1: Time sequences of an OE7 (d = 8 nm) vesicle fol- 
lowing electroporation. The first reassembled vesicle appears 
at t 26 sec, indicative of a large edge energy compared to 
the bending energy. The arrows indicate a particular disk 
reassembling. Scale bar is 10 /ion. 

For large pores (> 1 /itm) such as those in OE7 and 
OB16, we can visualize r(t) by optical microscopy. Al- 
though both OE7 and OB16 exhibit unstable pores, 
they both reassemble rather quickly (Fig. [j]). To our 
knowledge this is the first direct visualization of vesicle 
reassembly. Prior work has been focused in the subopti- 
cal regime and timescales for assembly have been of order 
minutes 0. This closure reflects the relative large mag- 
nitude of the edge energy (line tension) compared to the 
bending energy, characterized by a "vesiculation index" 
Vf ~ TRd j k c Jll| , where R<i is the radius of a disk and k c 
the bending modulus. By assuming an attempt frequency 
1/t z based on viscous drag in the solution |pl| , we can 

calculate for OE7, V) = 2 1 - ^/k b T ln(r Q /r z )/87rfc c « 

1.92. This reflects the rapid transition of OE7 into vesi- 
cles (i.e., small reassembly time r ) from what are pre- 
sumably disk-like or cylindrical micelles. 

Membranes of OB16 reassemble over a few minutes, 
indicating a smaller value of Vf compared to OE7. Vi- 
sualization of the OB 16 process was possible only with 
fluorescently labeled membranes, as the fragments drift 
substantially before reassembling. Initially, membrane 
is clearly lost as the pores grow to several microns in 
diameter (Fig. ||), again suggestive of a small line ten- 
sion relative to the surface tension. The hydrophobic 
dye indicates that OB16 has an increased preference 
to form cylindrical or disk-like aggregates upon loss of 
membrane integrity (Fig. ||). This is consistent with ob- 
servations that morphological boundaries (e.g., lamellar- 
to-cylindrical) shift toward lower hydrophilic block frac- 
tion with increasing M n [[l?], as well as with the longer 
reassembly time. OE7 and OB16 have characteristic 
reassembly times r a that are 20 and 100 sec, respec- 
tively. From scaling arguments pa], one can show that 



Figure 2: Time sequences of OB16 (d = 11 nm) following 
poration are visualized with fluorescently labeled membranes. 
Vesicle reassembles at t ta 110 sec. Note disk-like intermedi- 
ates with rough edges. Scale bar is 10 fim. 



T a ~ exp(c? 2 ), indicating that even if a larger copolymer 
such as OB19 were able to reassemble, the timescale 
would be several hours. Previous investigations of egg 
lecithin vesicle reassembly could not be visualized due 
to the rapid transition from disks to vesicles ]^8|. Only 
with the addition of an "edge-active" agent could disks be 
made sufficiently metastable for observation via electron 
microscopy. These findings are not surprising given that 
using the scaling relation above and taking d — 3 nm for 
a lecithin bilayer would correspond to r a of only a few 
seconds. 

By extending the analysis of hole growth in thin poly- 
mer films p9| , several groups have developed a frame- 
work for describing pore growth in cells and artificial 
vesicles There is an initial, viscous dissipation 

within the membrane which at a later time crosses over 
to a regime dominated by dissipation in the surrounding 
fluid. For OB16 (Fig. ||a), the only system we can unam- 
biguously observe, we find that the pore growth process 
is almost entirely described by viscous losses within the 
membrane. The time constant for the initial regime n 
(0.67 ± 0.05 sec) allows us to determine the membrane 
viscosity f] m = riS/d w 10 6 Pa sec. The resulting viscos- 
ity is three orders of magnitude larger than typical lipid 
membrane viscosities calculated by pore growth analysis 
(j] m « 10 3 Pa sec 0]), consistent with findings by other 
groups using different techniques ^0|. In the limit where 
the dissipation is primarily due to the surrounding fluid 
viscosity rji, pores will grow at a constant velocity | pl[ . 
From Fig. ||b, this is clearly not the case. We can also 
calculate the crossover radius r c to the inertial regime as 
Vmd/rji ~ 10 meters - indicating that this second regime 
is never reached for our micron-size vesicles. 

The small pores (< 1 fim) in OB 18 and OB 19 re- 




Figure 3: (a) Phase constrast imaging and (b) corresponding 
plot of pore radius r growth with time for OB16. The line 
is a theoretical fit that is dominated by viscous losses within 
the membrane. The resulting time constant n is 0.67 ± 0.05 
sec (N = 6). The error in measuring r is ±0.5 /im, limited by 
optical resolution. Scale bar is 10 /im 



Figure 4: (a) Loss of phase contrast in OB19 (d = 21 nm) 
due to escape of sucrose, (b) Corresponding normalized in- 
tensity of encapsulated sucrose versus time. Solid line is a fit 
to diffusion through a porous membrane, giving an effective 
permeability lo — 0.10 ± 0.03 ^tin/sec (N = 3). Scale bar is 
10 /xm. 



quire a slightly different analysis. In this case, the in- 
tensity of encapsulant I(t) is monitored, as pores cannot 
be directly observed (Fig. |^a). For a vesicle of size R v , 
there is an initial "leak-out" regime, where solute loss is 
driven by Laplace pressure AP = ZE/R V , lasting sev- 
eral milliseconds for nanometer-size pores. We cannot 
observe this regime, but the tension will relax on this 
timescale Tfc ~ rjiRl/rT,, and is followed by diffusion of 
solute through pores By fitting our data (Fig. fy>) to 
a simple exponential decay model of transport through 
porous membranes J2^|, we can obtain a time constant 
Td = -Rt>/3w, where u> is the effective membrane perme- 
ability. There is an initial time lag in the intensity decay 
tiag ~ 10 sec, which we tentatively attribute to hindered 
release of solute due to the PEO chains lining the nascent 



pores. The time constant increases with membrane thick- 
ness (rd = 13 ± 4 sec for OB18 and Td = 17 ± 6 sec 
for OB19), reflecting either slower dynamics within the 
membrane and/or smaller pore sizes. The net result is an 
effective (sucrose) permeability oj = 0.22 ± 0.15 /xm/sec 
for OB18 and w = 0.10 ± 0.03 /jm/sec for OB19. 

To estimate the size of sub-micron pores, we encapsu- 
lated fluorescent dextrans of various molecular weights 
and subjected these vesicles to electroporation (Fig. ||). 
Sucrose, which can be viewed under phase contrast and 
has a size ~ 0.6 nm, always escapes the porated mem- 
brane. The FITC-labeled dextran of M n w 4.4 kg/mol 
and a TRITC-labeled dextran of M„ w 160 kg/mol were 
retained over approximately 10 minutes. The sizes of 
these dextrans can be estimated by scaling relations |2^| 
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Figure 5: OB18 electroporation and encapsulant retention, (a) t < 0, sucrose appears dark under phase contrast microscopy, 
(b) At t = 1.2 min, phase contrast indicates that sucrose has escaped, while (c) fluorescence shows that TRITC-dextran (160 
kg/mol) is still retained. Scale bar is 25 /im. 



to give hydrodynamic radii of Rh sw 2 nm and Rh ~ 5 nm 
for the FITC and TRITC-dextrans, respectively. Thus 
although the pore size is only a few nanometers, given 
the dextran polydispersity (« 1.5), more precise mea- 
surements on the pore size will be published elsewhere 
& 

In summary, we have observed for the first time vesicle 
formation from disk-like fragments, allowing us to obtain 
scaling relations for reassembly time. In certain cases, 
the slow dynamics of pore growth can be used to ob- 
tain the membrane viscosity (e.g., OB16 i] m — 10 6 Pa 
sec), qualitatively consistent with observations by other 
groups ^0|. Sucrose leakage from thicker membranes 
(OB18 and OB19) gives effective membrane permeabil- 
ities of vesicles with nanometer-scale pores. Control of 
the size and number of these pores by methods such as 
crosslinking |2^| invites speculation on possible applica- 
tions such as controlled drug delivery as well as for experi- 
mental tests of predictions for (bio) polymer translocation 
through pores [^5[ p6| . The results illustrate the rich di- 
versity of behavior possible by changing a critical length 
scale d, with other chemistries potentially providing even 
more distinctive effects. 
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